Impurity Effects in Sign-Reversing Fully Gapped Superconductors: 
Analysis of FeAs Superconductors 
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To understand the impurity effect on Tc in FeAs superconductors, we analyze a simple two-band 
BCS model with a repulsive interband interaction. The realized fully gapped superconducting state 
with sign reversal, which is predicted by spin fluctuation theories in this compound, is suppressed 
by impurities due to the interband hopping of Cooper pairs, if the interband impurity scattering 
I' is equal to the intraband one /. When |/'//| < 1, in high contrast, Tc is almost unchanged by 
strong impurity scattering since interband scattering is almost prohibited by the multiple scattering 
effect. Since \I' /I\ ~ 0.5 is expected, the robustness of superconductivity against impurities in FeAs 
superconductors is naturally understood in terms of the sign-reversing fully gapped state. 
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After the discovery of the superconductor 
La(0,F)FeAs (Tc = 26 K) [1], high-Tc supercon- 
ductors with FeAs layers have been studied intensively 
[2-5]. The ground state of the undoped compound is 
the spin density wave (SDW) state, where the ordered 
magnetic moment is ^ 0.3 /ib and the ordering vector is 
Q « (7i',0) or (0,7r) [6-9]. The superconducting state is 
realized next to the SDW state by carrier doping [9, 10] 
According to NMR study, the singlet superconducting 
state is realized in FeAs [11-13]. In the first-principle 
band calculations [14, 15], the Fermi surfaces in FeAs 
are composed of two hole-like Fermi pockets around 
the r = (0, 0) point (FSl and FS2 in Fig. 3) and two 
electron-like Fermi pockets around M= (tt, 0), (0, tt) 
points (FS3 and FS4 in Fig. 3). 

Theoretically, there are several possible pairing states: 
According to the random phase approximation (RPA) 
based on a realistic five-orbital tight-binding model [16], 
the nesting between the hole and electron pockets gives 
rise to the strong antiferromagnetic (AF) fluctuations 
with Q « (tt, 0), which is consistent with experimental 
results. Then, a fully gapped s-wave state with sign re- 
versal is expected to emerge since AF fluctuations works 
as the repulsion interaction between hole and electron 
pockets [16-21]. Although AF fluctuations due to the 
nesting between two electron pockets [q ~ (tt, 7r/2)] can 
induce the d^2_y2--wave state with line nodes on the hole 
pockets, the obtained Tc is rather low [16, 18-20]. On 
the other hand, the conventional s-wave state without 
sign reversal will be realized if the charge fluctuations or 
electron-phonon interactions are strong. 

Experimentally, a fully gapped superconducting state 
has been determined by recent penetration depth mea- 
surement [22], angle-resolved photoemission spectroscopy 
(ARPES) [23-25], and specific heat measurement [26]. 
The resonance peak observed by inelastic neutron mea- 
surement below Tc supports the sign-reversing super- 
conducting gap mediated by AF fluctuations [27, 28]. 
Anomalous transport phenomena (such as Hall coeffi- 



cient and Nernst signal) in the normal state in FeAs, 
which are similar to those observed in high-Tc cuprates 
and CeAflns (M=Co, Rh, Ir) [29], also indicate the ex- 
istence of strong AF fiuctuations [10, 12, 30, 31]. On 
the other hand, impurity effect on Tc due to Co, Ni, 
or Zn substitution for Fe sites is very small or absent 
[11, 12, 32-34], which is decisive for determining the pair- 
ing symmetry. For example, this result clearly rules out 
the possibility of line-node superconductivity. This re- 
sult may also eliminate the theoretically predicted sign- 
reversing s-wave states, since the Cooper pair is de- 
stroyed by the interband scattering induced by impu- 
rities; (k t,-k Dbandi (k' T,-k' i)band2- In this 
manner, a study of the impurity effect on FeAs super- 
conductors offers us significant information on the pair- 
ing symmetry, and therefore reliable theoretical analyses 
are highly required. 

In this letter, we study a simple two-band BCS model 
to investigate the impurity effect on FeAs superconduc- 
tors. Therein, the sign-reversing pairing state is realized 
if we introduce the interband repulsive interaction to de- 
scribe the effective interaction due to AF fluctuations. 
When the interband impurity scattering potential /' is 
equal to the intraband one /, the reduction in Tc per im- 
purity concentration ni,„p, — ATc/nimp, is prominent as 
in non-s-wave superconductors. However, x = \I' /I\ is 
smaller than 1 in this compound since hole and electron 
pockets are not composed of the same d-orbitals. In this 
case, — ATc/riimp becomes very small; in particular, it 
approaches zero in the unitary regime (//Wband ^ 1)- 
Therefore, the experimental absence of the impurity ef- 
fect on Tc in FeAs is well understood in terms of the 
sign-reversing s-wave state proposed in refs. [16] and 
[17]. 

Here, we analyze the Eliashberg gap equation for the 
two-band BCS model. To concentrate on studying the 
impurity effect, we neglect the mass-enhancement factor 
and the quasiparticle damping due to electron-electron 
interaction for simplicity, both of which are given by the 
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normal self-energy. In the absence of impurities, the lin- 
earized gap equation at Tc is given by [35]: 



Ai{pn) 



(1) 



j=a,l3 



where PmPm are fermion Matsubara frequencies, and 
i,j (= a,/3) are band indices. Ni is the density of states 
at the Fermi level, and gij is the effective interaction 
between band i and band j. = ^("^c ~ bm|), 

where lOc is the characteristic energy scale of the effec- 
tive interaction. Ai(p„) is the gap function, which is 
independent of p„ for |p„| < lOc in the absence of im- 
purities. Hereafter, we assume gap = g > to real- 
ize the sign-reversing s-wave gap (A^A^ < 0), and put 
gaa = 3/3/3 = for simplicity since they will be much 
smaller than g in FeAs. After the standard analysis [35], 
the transition temperature without impurities is obtained 
as = 1.13wcexp(-l/.g^iV„iV/3). 
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FIG. 1: (a) T-matrix in the normal state, (b) Impurity- 
induced normal self-energy (iJE^) and the linealized anoma- 
lous self-energy {ST/}^) in the T-matrix approximation. 

Hereafter, we study the nonmagnetic impurity effect 
in the two-band BCS model with g^p > 0. In similar 
models, the effect of interband impurity scattering in the 
NMR relaxation ratio has been analyzed using the Born 
approximation [36] and T-matrix approximation for / = 
I' = oo [37]. Very recently, Bang et al. have reported 
a sizable reduction in Tc by a strong impurity potential 
with 1 = 1' [38]. In contrast to their result, we will show 
below that Tc is almost unchanged when /'// < 1. 

Hereafter, we use the T-matrix approximation, which 
gives the exact result for riimp <C 1 for any strength 
of We can assume that > without los- 

ing generality. Using the local normal Green function 
9t{Pn) = -iirN^Sn {sn = sgn(p„)) [35], the T-matrix in 
the normal state is given by 



E 

l=a,p 



higiiPn)Tlj{pn), 



(2) 



where = ISij + — 5ij). Its diagrammatic expres- 
sion is shown in Fig. 1 (a). Except at x = 1, the solution 
of eq. (2) behaves as Taa ^ —iSn/T^Na + 0{I^^) and 
Tap ~ in the unitary limit. This fact means that 

the the superconducting state is unaffected by a strong 



impurity potential since the interband scattering is pro- 
hibited as a result of multiple scattering. 

Using the T-matrix, the impurity-induced normal and 
anomalous self-energies are respectively given as 

(5E"(p„) = riimpTiibn), (3) 

5T.1{pn) = riimp ^ Tu{pn)fl{Pn)Tu{-pn), (4) 
l=a,P 

where fi{pn) = TTNiA.i{pn) /\pn\ is the local anomalous 
Green function at Tc. Their expressions are shown in 
Fig. 1 (b). Then, the gap equation at Tc is given by [35] 

Zi{pn)Ai{pn) = -TtTc ^ gijNj^ Aj{p.m)/\Pni\ 
j=a,P m 

+'5S-bn), (5) 

where Z,{p„) = 1 - (SJ^UPn) - <5Ef (-p„))/2^p„ = 1 + 
"fi/\pn\, and 7i = -ImEf (p„) • s„ (> 0) is the quasipar- 
ticle damping rate due to impurity scattering. 

Now, we analyze the impurity effect on Tc in the case 
of Na = Np = N a.s the first step. In this case, the 
relationships ja = ip = 1, Aa{pn) = -Ap{pn) = A(p„), 
and (JSJ^ = ~<^^/3 = satisfied. They are obtained 

as 

7 - n,^p7:N[{l' + I'') + n^N\l' - I'')^]/A, (6) 
S^'^iPr.) = ni^pA(p„)7r7V[(/2 - 

+n'N\l'-I''y]/\p,,\A, (7) 



where A is defined as 



A = l + 27r2iV2(/2 + /'^) + n^N\l^ - l") 



(8) 



We will show below that the interband impurity scatter- 
ing /' is renormalized by l/vC4. Using eqs. (5)-(8), the 
frequency dependence of the gap function is obtained as 



A(p„) = C[Z(p„)~'5S''(p„)/A(p„)]"^ 



C 



l-|-2nu„p7riV/'71p„l^ 



(9) 



where C = t:N gTc^[^ A{pm) /\Pm\ is a constant inde- 
pendent of pn. By inserting eq. (9) into the definition of 
C, we obtain the following equation for Tci 



1 = ^iV^Tc^' [lp™l + 27rni,„pArj''M 
gN In- 



27rTc 




(10) 



where ^p{x) is the digamma function. The equation for 
Tc° is given by dropping the term UimpNI' jT^A in eq. 
(10). It is noteworthy that this term also vanishes when 
X ^ \ and / — > oo, even if riimp > 0. This fact means 
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that superconducting state is unaffected by impurities in 
the unitary limit. 

Using eq. (10) and the relation ip'{l/2) = 7r^/2, we ob- 
tain the relationship ln{T^ /Tc) = nin,pTT'^Nr'^/2T°A for 
<C 1. As a result, the reduction in Tc per impurity 



' ''imp 

concentration for 



^ 1 is obtained as 



2A 



(11) 



The physical meaning of the right-hand side of eq. (11) 
is the rate of pair breaking, which is given by the ampli- 
tude of interband scattering for Cooper pairs: ITq,^!^. Its 
I dependence is shown in Fig. 2 (a). When x = 1, eq. 
(11) increases in proportion to in the Born regime 
{ttIN <^ 1), and it approaches 1/8A^ in the unitary 
regime {ttIN ^1). In the latter case, the superconduc- 



tivity in FeAs will be destroyed only at n 



imp 



0.02 since the average between the electron and hole den- 
sity of states per Fe atom is 0.66 eV~^ [14]. When x ^ 1, 
in contrast, eq. (11) ^ x"^ /2tt'^N^P{1 - x^f ^ in the 
unitary limit. In this case, pair breaking is almost absent 
and Tc is unchanged. 
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FIG. 2: (a) -ATc/nimp and (b) -Arc/7 as functions of /. 
l/8iV ~ 2500 K in FeAs. When x < 0.8, -ATc becomes very 
small for / > 3. 



According to eqs. (11) and (6), we obtain as 



AT, 



7 2(l + a;2) + 27r2iV2/2(i_a;2) 
In the Born limit 



(12) 



-Ar,/7 



/2(1 + a;2), which 



diminishes slowly as x decreases from unity. In the 
unitary limit, it is strongly suppressed as — AT^/t ~ 
[- ATc /7] Born X {nNI)-^ for a: ^ 1. Note that 7 
is related to residual resistivity as po = 2m'y / e^n, where 
n is the carrier density. Figure 2 (b) shows the / de- 
pendence of eq. (12): For x = 0.5, -ATc/j = 0.0077 
at IN ~ 3, which is two orders of magnitude smaller 
than 7r/4. Such a small impurity effect will be difficult 
to observe experimentally. 

Next, we discuss the case of Na 7^ N/s, where the rela- 
tion Aa = — A^ is not satisfied. Even in this case, we can 
obtain — ATc/nimp by solving eq. (5) analytically. After 
a long calculation, the obtained result for riimp <C 1 is 



3(A^„+iV;3)-2^ArA^]/'' 



8A 



(13) 



1 + 7r^P{N^ + 



2NaNf3TT^r^ 



/'■^)2, which is proportional to — x"^) 



2\2 



ATe ^ 

where A = 
NINIt:\P - 

in the unitary regime. Therefore, eq. (13) approaches 
zero in the case of a; 7^ 1 in the unitary regime. Thus, 
the obtained results for Na — Np given in Fig. 2 are 
qualitatively unchanged even for Na 7^ A'^- 
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FIG. 3: Fermi surfaces (FSs) in FeAs superconductors in the 
unfolded Brillouin zone; see ref. [16]. 

On the basis of the above results, we discuss the im- 
purity effect on Tc in real FeAs superconductors. Figure 
3 shows the Fermi surfaces (FS1-FS4) of this compound 
in the unfolded Brillouin zone [16, 17]: FS1,2 are mainly 
composed of dxzidyz orbitals of Fe, whereas FS3,4 are 
composed of dxzi dyz and dx^-y2 orbitals, according to the 
five-orbital model in ref. [16]. More precisely, FS3a,4a 
(3b,4b) are composed of d^z , dyz orbitals {dx'^-y^ orbital) . 
Mazin et al. had proposed that superconducting gap 
functions in FS1,2 (Ai,2) and those in FS3,4 (A3^4) are 
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different in sign. In this case, bands a and /? in the 
present study correspond to FS1,2 and FS3,4 in FeAs, 
respectively. 

Here, wc consider an impurity rf-atom (such as Co, Ni. 
or Zn) placed at an Fe site. In the d orbital representa- 
tion, the local impurity potential will be diagonal with 
respect to the d orbital d{d!y. {I)d.d' = ISd.d'- In the 
band-diagonal representation, it is transformed into 

« /(X^Od,i(k)Od,i(k'))^L,k'e/3, (14) 

d 

where i,j (= a,f3) represent the band indices. Od.i{^) = 
(d;k|z;k) is the transformation orthogonal (or imitary) 
matrix between the orbital representation (orbital d) and 
the band-diagonal representation (band i). Therefore, 
laa, Ipp ~ I whereas |/a/3l should be smaller than |/| 
when bands a and /3 are composed of different d orbitals. 
Since FS3b,4b are composed of d^2_y2 orbitals of Fe, the 
relation x = \I' /I\ ~ 0.5 is realized in FeAs. 

On the other hand, according to the RPA analysis by 
Kuroki et al. [16], A3. 4 has line nodes near (tt, 0) — (0, tt) 
line owing to spin fluctuations with q « (tt, 7r/2); therein, 
the sign of for FS3a,4a (3b, 4b) and that of Ai,2 
are equal (different). Then, bands a and /3 in the 
present study correspond to [FSl,2-|-3a,4a] and [FS3b,4b] 
in FeAs, respectively. Since bands a and /3 are composed 
of different d orbitals, the relation .t <C 1 is expected. 
Note that the obtained s-wave state is fully gapped since 
the line nodes on FS3,4, which arc not protected by sym- 
metry, arc masked by Kniall intcrband pairing [16]. We 
also note that the sign change in A3^4 does not occur and 
Mazin's type sign-reversing state is realized in the case 
ofUfiiU' [16]. 

As a result, in both unconventional s-wave states pro- 
posed in refs. [16] and [17], -ATc/7™p « in the uni- 
tary regime. In high-Tc cuprates, / due to Zn impurity is 
about 10 eV [31]. If we expect / ~ 10 eV and N^,/} ~ 1 
eV~^ in FeAs, the unitary regime is actually achieved. 

Although we have put rjaa = fjpp = above, they will 
take positive (negative) values owing to electron-electron 
correlation (electron-phonon interaction) in FeAs super- 
conductors. Although these diagonal interactions modify 
T°, we can show that both eqs. (11) and (12) are un- 
changed even if gaa,9i30 0- Thus, the obtained results 
in the present study will be valid for general multiband 
fully gapped superconductors with sign reversal. 

In summary, we analyzed the impurity effect on Tc in 
a sign-reversing s-wave BCS model. Except at I' = I , Tc 
is almost unaffected by impurities in the unitary limit, 
since the interband elements of the T-matrix vanish due 
to multiple scattering. Thus, the robustness of super- 
conductivity against Co, Ni, or Zn impurities in this 
compound is naturally explained in terms of the sign- 
reversing s-wave superconductivity [16-21]. On the other 
hand, a weak impurity scattering causes pair breaking. 



Therefore, the finite density of states may be induced in 
the superconducting state by weak impurities or disorder. 
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